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We introduce a correlated static model and investigate a percolation transition. The model is a 
modification of the static model and is characterized by assortative degree-degree correlation. As one 
varies the edge density, the network undergoes a percolation transition. The percolation transition 
is characterized by a weak singular behavior of the mean cluster size and power-law scalings of the 
percolation order parameter and the cluster size distribution in the entire non-percolating phase. 
These results suggest that the assortative degree-degree correlation generates a global structural 
correlation which is relevant to the percolation critical phenomena of complex networks. 
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I. INTRODUCTION 

Many systems in nature have a complex network struc- 
ture. For example, the Internet is a wired network of 
computers and routers, the World Wide Web is a network 
of web pages hyperlinked to each other, a protein inter- 
action network is a network of interacting proteins in an 
organism, and a social network is a network of individuals 
who are linked through a certain relationship. Such com- 
plex networks have a highly inhomogeneous structure. In 
order to characterize and understand the structure and 
dynamics, extensive research has been performed for the 
last decade 0, fl & |, i, 1, . 

A degree distribution p(k) for the probability of a node 
having degree k is a quantity of primary importance. It is 
one of the essential characteristics that influences struc- 
tural properties, dynamical behaviors, and collective phe- 
nomena of complex networks. However, the degree is a 
property of each individual vertex. Hence, the degree 
distribution by itself cannot describe correlation among 
different vertices. It is attracting growing interest, since 
many real-world networks display a certain level of cor- 
relation, which has a significant impact on network prop- 
erties H 1 M EH • 

The degree-degree (DD) correlation refers to the cor- 
relation between degrees of neighboring vertices [1] . The 
whole information on the correlation is contained in the 
degree correlation function p{kl , k), the probability of an 
edge linking nodes of degree k and k', and the con- 
ditional probability p(k' \ k) = p(fc', fc)/(^ fc „ fc)), 
the probability of a node among neighbors of degree-fc 
nodes having degree k' . The overall feature is conve- 
niently characterized by the assortativity and the mean 
neighbor degree function. The assortativity coefficient r 
is defined as the normalized Pearson correlation coeffi- 
cient between degrees of neighboring vertices @ , and the 
mean neighbor degree function is defined as K^Nik) = 
J2k' k'p(k'\k) A network with a positive (neg- 
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ative) value of r or an increasing (decreasing) function 
ifjvjv(^) has a positive (negative) correlation and is called 
assortative (disassortative). A network with r = or a 
constant function KNN{k) is called uncorrelated or neu- 
tral. 

Recent studies have shown that the structural corre- 
lation, as well as the degree distribution, plays an im- 
portant role in percolation critical phenomena (ill . 
In Ref. , the percolation transition was studied in the 
exponential random graph (ERG) model. With the ERG 
model one can simulate assortative, neutral, or disassor- 
tative networks with the same degree distribution. The 
numerical study reveals that the disassortative and the 
neutral ERG models display the percolation transition in 
the same universality class. On the other hand, the as- 
sortative ERG model displays the percolation transition 
in a distinct universality class. The percolation transi- 
tion in the assortative ERG model is similar in nature to 
that in the growing-network models (Til . [3113,111 which 
are also assortative. It is noteworthy that the percola- 
tion transition in the assortative ERG model cannot be 
explained only with the DD correlation A possible 
scenario is that the assortative DD correlation may give 
rise to a global structural correlation that is relevant to 
the percolation critical phenomena. 

In this work, we introduce a model for assortative net- 
works and investigate the percolation transition of the 
model. The purpose is to confirm whether the assorta- 
tive DD correlation is relevant to the percolation critical 
phenomena. Another purpose is to present an efficient 
stochastic model of correlated networks with degree dis- 
tribution ranging from the Poisson distribution to scale- 
free power-law distributions. The ERG model introduced 
in Ref. [lj] is based on a Monte Carlo method which takes 
a longer time to generate large-size networks. We present 
a model with which one can generate a correlated network 
fast and easily. This model can also be used for further 
studies of critical phenomena other than the percolation. 

In Sec. HH we introduce a model for networks with an 
assortative DD correlation. This is a modification of the 
so-called static model [2(| which is a model for uncor- 
related scale-free networks. Our model will be referred 
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to as the correlated static model. Basic properties of the 
model are also presented. In Sec. lIIIi we investigate bond 
percolation transitions in the correlated static model. We 
summarize and conclude the paper in Sec. HVl 



II. CORRELATED STATIC MODEL 

The static model is an efficient model for uncorre- 
cted networks [20] . A static- model network with N ver- 
tices and K edges is constructed as follows: (i) Each 
vertex i (i = 1, • • ■ , N) is assigned to a weight Wi = 

i _A V(Sjli J ( u ) according to the probability {wi}, 
two vertices are chosen at random. If there is no edge be- 
tween them, they are linked with an edge. The procedure 
(ii) is repeated until one has K edges in total. A resulting 
network is scale- free with a power-law degree distribution 
p(k) ~ k~ x with the degree exponent A = 1 + 1/ fj,. When 
fj, = 0, all vertices are chosen with equal probability and 
the model reduces to the Erdos-Renyi random network 
with the Poisson degree distribution. The static-model 
network is uncorrelated for \i < 1/2 or A > 3 [21]. The 
percolation transition in this model has been thoroughly 
studied [22T ]. 

We modify the static model to incorporate the assor- 
tative DD correlation. A network with A vertices and 
K edges in the correlated static model is constructed as 
follows: (i) Each vertex i (i = 1, • • • , A) is assigned to 

a weight Wi = £ _/ V(X)j=i •? '~^)-- ( n ) a vertex * is cho- 
sen with probability W{ at random, and then another 
vertex j is chosen randomly among all vertices with the 
same degree as i. If there is no edge between % and j, 
an edge connecting them is added. The procedure (ii) 
is repeated until there are K edges in total. Note that 
the static model and the correlated static model differ in 
the procedure (ii). While edges are added between ver- 
tices chosen independently in the former, they are added 
between vertices of the same degree in the latter. This 
generates a positive DD correlation. 

Consider the degree distribution of the correlated 
static model. Let pi(k,t) be the probability that a 
vertex i has degree k at time step t. Since an edge 
is added at each time step, the total number of edges 
K is equal to t. The degree distribution is given by 
p{k,t) — jj^iPi{k,t). For the sake of convenience, we 
introduce qk(t) = Y^i w iPi(^i^)- Then, in the large-A 
limit, the time evolution of pi(k,t) is governed by 



Pi(k,t+ 1) -pi(k,t) 



gfc-iW 
Np k -i(t) 

Qk(t) 
Np(k,t) 



Pi(k-l,t) 
Pi(k,t), 



with the initial condition Pi(k,t = 0) = S(k,0). The 
terms in the square brackets become 2uii for the static 
model. Summing up both sides over all i, one finds 
the evolution equation of the degree distribution func- 




FIG. 1: (a) Degree distributions of the static model (solid 
line) with N — 2 9 x 10 3 and the correlated static model with 

N = 10 3 (■ ■ ■), 2 3 x 10 3 (- - -), 2 6 x 10 3 ( ), and 2 9 x 10 3 (- 

•— ) with n = 0.5 and link density 5 = K/N = 2.0. (b) 
Scaling of the maximum degree at fi = 0.2, 0.5, and 0.8. 
Open symbols are for the static model and closed symbols 
are for the correlated static model. Solid lines have slope 
l/\ = (i/(l + fi) and dotted lines have slope 1/(A — 1) = fi. 



tion p(k, t). This is given by 

p(k, t + l)= p{k, t) + A [q k -!(t) - q k {t)} . (1) 

Note that the correlated static model has the same evo- 
lution equation as the static model. Hence we conclude 
that the correlated static model has the same degree dis- 
tribution as the static model. From the properties of the 
static model [20j . follows the power-law degree distribu- 
tion 



p(k) 

with the degree exponent 



A = 1 + l/(i 



(2) 



(3) 



for [i > and the Poisson distribution at /i = 0. 

We present the degree distributions of the correlated 
static model and the static model with /i = 0.5 in 
Figure [TJa). This confirms the theoretical prediction 
that both models have the power-law degree distribu- 
tion with the same degree exponent A. However, we 
find that the maximum cutoff degree k max scales dif- 
ferently with the network size N . The maximum de- 
gree k m ax is known to scale as k max ~ A 1 /^ 1 ) in the 
static model [23]. This can be derived from the condi- 
tion that the number of nodes with k > k max should 
be of the order of unity. This condition demands that 

■CL*P( fc ) dfc ~ fZaa, k ~ X<lk = <D ( 1 / N )' WhiGh l0adS t0 

k m ax ~ A 1 /^ -1 ). On the other hand, as one can see 
in Figure DIb), k ma x in the correlated static model is 
much smaller than that in the static model. The fol- 
lowing argument explains the scaling behavior of k max . 
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FIG. 2: (a) Assortativity coefficient of the static model (open 
symbols) and the correlated static model (closed symbols), 
(b) Knn vs. k in the static model (lines) and the correlated 
static model (symbols) with N = 512000 vertices. 



In the correlated static model, one can add a link to a 
vertex only if there exists another node with the same de- 
gree. So, if a node has the maximum degree k max , there 
should be another node having the same degree k max . 
This leads to the constraint p(fc m ax) = 0(1/N) instead 
of p(k)dk = 0{1/N). This gives 



(4) 



The numerical data in FigureQJb) are in good agreement 
with the scaling behavior in Eq. ((3]). 

In order to examine the DD correlation, we measure 
the assortativity coefficient r [8|. This is defined as 



(fc e ,lfc e ,2>e " ((fc e ,l + fc e , 2 )/2} 2 



((*; 



e.l 



fce 2 2 )/2} £ 



((fce,l + k e , 2 )/2)l 



(5) 



where ((-)) e represents the average over all edges, and 
k £y i and fc ej 2 denote the degrees of two vertices connected 
with an edge e. It is measured and presented in Fig- 
ure HJa) . While the assortativity coefficients vanish as 
TV increases in the static model, they converge to a finite 
value as N increases in the correlated static model. Pos- 
itive correlation is also observed in the mean neighbor 
degree K^Nik) presented in Figure HJb). In the cor- 
related static model Kj^iq(k) ~ k, while K^N(k) is a 
constant function in the static model. This shows that 
the correlated static model has the desired property. 



III. PERCOLATION TRANSITION 

As a static-model network does [22j], so also a cor- 
related model network undergoes a percolation transi- 
tion as one increases the edge density 6 = K/N. To a 
given value of 5, vertices are decomposed into disjoint 
sets (called clusters) in such a way that all vertices in a 
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FIG. 3: (a) S vs. 8 for the correlated static model (a) and the 
static model (b) at /i = (solid lines) and 0.2 (dotted lines). 
Network sizes are N = 10 5 , 4 x 10 5 , 16 x 10 5 , and 64 x 10 5 . 
The larger N is, the higher the peak is. (c) Peak height Smax 
vs. N. Closed symbols are for the correlated static model, 
and open symbols are for the static model. The solid line has 
slope 1/3. 



cluster are mutually connected, while those in different 
clusters are not. The size s of a cluster is defined as the 
number of vertices in it. When there is no edge (5 = 0), 
each vertex belongs to a cluster of size s = 1. As edges 
are added, small clusters merge into large ones. A cluster 
configuration can be characterized by the cluster size dis- 
tribution n s , which is defined as the ratio of the number 
of clusters of size s to N. 

The order parameter for the percolation transition is 
given by 



p = 

1 no — 



N 



(6) 



where s max is the size of the largest cluster. A network 
is said to be in a percolating phase if is finite in the 
— > oo limit. In that case, the largest cluster is called 
the infinite or giant cluster. If vanishes in the A^ — > oo 
limit, a network is said to be in a non-percolating phase. 
Another useful quantity is the mean cluster size S, given 
by 



5' 



s 2 n s 



(7) 



where the summation is over all clusters except for the 
largest cluster. 

In the static model, the percolation transition is char- 
acterized by power-law scalings, ~ (6 — 6 C )P and 
S - \5-5c\-i for A > 3, where S c = (1 - 2/*)/ (2(1 - /i) 2 ) 
is the percolation threshold. The critical exponents have 
the values (3 = 1 and 7=1 for A > 4, and j3 = 1/(A — 3) 
and 7 = 1 for 3 < A < 4 3 . When A < 3, the networks 
are always in the percolating phase. 

We perform numerical studies of the percolation tran- 
sition in the correlated static model and compare the 
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FIG. 4: Plot of the effective exponent a e (N,S) in the corre- 
lated static model with fi — 0.0 (a) and /U = 0.2 (b), and in 
the static model with fx — 0.0 (c) and fx — 0.2 (d). Circle, 
square, and diamond symbols are used for N — (1, 4, 16) x 10 , 
respectively. 

results with those of the static model. A striking differ- 
ence is observed in the behavior of the mean cluster size 
S. In Figure 02 we compare S obtained for the static and 
correlated static models at fj, = and 0.2. In the static- 
model cases, there are peaks in the plot of S near per- 
colation thresholds, which sharpen as N increases. The- 
oretically, the peak heights should grow algebraically as 
Smax ~ N 1 ^ with the finite-size-scaling exponent v = 3 
for A > 4 and v = (A - 1)/(A - 3) for 3 < A < 4 0. 
Our data are consistent with the power-law scaling for 
the static model. 

On the other hand, in the correlated-static-model 
cases, the peaks are not as sharp as in the static model. 
The peak heights seem to converge to a finite value or 
at most scale logarithmically with N. The weak singular 
behavior of S suggests that the percolation transition in 
the correlated static model does not belong to the same 
universality as that in the uncorrelated static model. 

In order to characterize the phase transition, we study 
finite-size-scaling behaviors of the order parameter P^ 
near percolation thresholds. At criticality, we expect that 
Poo scales algebraically as 

P^N) ~ N- a . (8) 

The asymptotic value of the scaling exponent a is 
obtained from the analysis of the effective exponent 
a e (N,5) = -ln(P 00 {mN,S)/P 00 (N,S))/hxm with con- 
stant m = 4. If Pqo follows the power law, then the 
effective exponent will converge to the scaling exponent 
in the N — > oo limit. Otherwise, it will converge to a 
trivial value of or 1. 

In Figure QJ we present a plot of a e (N, S) for the static 
model and the correlated static model. We find a sig- 
nificant difference in finite-size-scaling behaviors of a e . 
First, consider the correlated static model shown in Fig- 
ure HJa) and (b). At large values of S, a e vanishes as N 
increases, which indicates that the system is in the perco- 




FIG. 5: Cluster size distribution, n s , in the non-percolating 
phase of the correlated static model with /i = 0.0 (a) and 
0.2 (b). 



fating phase with finite P*,. On the other hand, at small 
values of S, it converges to a nontrivial value a < 1 that 
varies continuously with 5. This suggests that the corre- 
lated static model is critical not only at the percolation 
threshold but also in the entire non-percolating phase. 
We estimate the percolation threshold as the point at 
which the power-law scaling sets in. The results are 
S c = 0.19(1) at n = 0.0 and S c = 0.15(1) at fi = 0.2. 
At the transition point, the scaling exponent takes the 
value 

a c = 0.5(1) (9) 

in both cases with \i = 0.0 and 0.2. 

The unusual scaling behavior of Poo in the correlated 
static model becomes clear when one compares it with 
that in the static model. Figure QJc) and (d) show that 
the static-model networks are critical only at the perco- 
lation thresholds, S c ~ 0.50 at fx = 0.0 and S c ~ 0.47 
at /i = 0.2. The scaling exponent at the critical point is 
given by a c ~ 0.33. These numerical results are in good 
agreement with the analytic results (22| . 

In order to investigate further the criticality in the non- 
percolating phase, we study the cluster size distribution, 
n s . Figure O shows n s at 6 < S c in the correlated static 
model with fj, = 0.0 and 0.2. We find that the cluster size 
distribution follows a power-law distribution n s ~ s~ T 
with varying exponent r. This result also suggests that 
the correlated static model is critical in the entire non- 
percolating phase. 

We also studied the percolation transition at larger 
values of /i (smaller value of A). In the static model, 
networks are always in the percolating phase (<5 C = 0) if 
fi > 1/2 or A < 3 [22|]. However, our numerical data show 
that the percolation threshold in the correlated static 
model already vanishes at y, = 1/4 (A = 5). This implies 
that the correlated static model network is more robust 
than the static model network. 
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IV. SUMMARY 

We have introduced a correlated static model and in- 
vestigated the nature of the percolation transition. Inter- 
estingly, the correlated static model is critical in the en- 
tire non-percolating phase below the percolation thresh- 
old. That is, the percolation order parameter scales alge- 
braically as Poo ~ N~ a and the cluster size distribution 
follows the power law as n s ~ s _r with varying exponents 
a and r. Such a type of percolation transition is also ob- 
served in the growing-network models and the assortative 
exponential random-graph model. These results support 



the claim that the assortative DD correlation may give 
rise to a structural correlation which is responsible for 
novel type of percolation transition. 
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